4D-XY quantum criticality in a doped Mott insulator 
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A new phenomenology is proposed for the superfluid density p s of strongly underdoped cuprate 
superconductors based on recent data for ultra-clean single crystals of YBa2Cu307_2,. We show that 
the puzzling departure from Uemura scaling and the decline of the slope as the T c — quantum 
critical point is approached can be understood in terms of the renormalization of quasiparticle 
effective charge by quantum fluctuations of the superconducting phase. We then employ (3+1)- 
dimensional XY model to calculate, within particular approximations, the renormalization of p s 
and its slope, explain the new phenomenology, and predict its eventual demise close to the QCP. 
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The manner in which superconductivity in high-T c 
cuprates gives way to Mott insulating behavior is a long- 
standing puzzle of fundamental importance. The anoma- 
lous behavior is revealed most strikingly in studies of the 
doping (x) and temperature (T) dependence of the su- 
perfluid density p s 1, 2]. As the doping is reduced, both 
p s and the critical temperature T c decline while the maxi- 
mum superconducting gap Ao increases. This dichotomy, 
along with the empirical Uemura relation Trjxp s [3] , led 
to suggestions by numerous authors 0, USE Hi that the 
superconducting transition in the underdoped cuprates 
is a phase- disordering transition and the pseudogap state 
above T c should be thought of as a "phase-disordered" 
d-wave superconductor. Ex per iments indeed show evi- 
dence for magn etic vortices [Tj|, strong fluctuation dia- 
magnetism ^l[ , and fermionic nodal quasiparticles in 
the pseudogap state of high- J), compounds. In addition, 
there is credible evidence that the superconducting tran- 
sition in many compounds is in the 3D-XY universality 
class [lj, [lj, [Tj, [l(| with a wide critical region, exactly 
as one would expect near a phase-disordering transition. 

The focus of the present Letter is to obtain a theo- 
retical understanding of the behavior of p s in strongly 
underdoped cuprates. Theoretical efforts to date have 
overwhelmingly addressed the phenomenology summa- 
rized by Lee and Wen 0. However, recent experiments 
performed in unprecedented proximity to the Mott in- 
sulator in ultraclean single crystals of YBa2Cu307_a; 

(YBCO) mEm 

and in high quality films |20j are qui- 
etly overturning the old paradigm: (i) Unlike optimally 
doped and weakly underdoped samples, the strongly un- 
derdoped data show no visible 3D-XY critical region in 
samples with T c < 25K; rather, the approach to T c is 
mean-field like, (ii) Measurements indicate a relation- 
ship between the critical temperature and the T = 
superfluid density of the form 



3d-xy thermal 



T c (xp], with 7-0.4-0.7, 



(1) 



a significant departure from the Uemura scaling. (Hi) 
The overall doping and temperature dependence can be 
parametrized as 



Ps (x,T)~ Ay 2 -By(k B T), 



(2) 




where y = T c (x)/T™ ax is a measure of the doping x, 
ymax _ ggj^ j g maximum critical temperature for 



FIG. 1: Behavior expected in the vicinity of the (3 + z)D- 
XY quantum critical point in the doping-temperature plane. 
Solid line T c (x) represents the superconducting phase transi- 
tion; dashed lines are crossovers. 



YBCO, A ~ 66meV, and B ~ 9.5. The demise of su- 
perconducting order described by Eq. (|2J signals a pro- 
found departure from predictions of RVB-type theories 
and earlier parametrizations of p s (x,T) [2I l2l|. In what 
follows, we demonstrate how this new phenomenology 
follows simply and elegantly from an effective theory de- 
scribing a phase-fluctuating d-wave superconductor. 

The absence of any visible 3D-XY fluctuation region 
in the data is surprising since one would naively expect 
phase fluctuation effects to become more pronounced in 
the underdoped region as the system approaches the 
Mott insulator. Upon closer inspection, however, one 
finds that the observed behavior is entirely consistent 
with the behavior of a system approaching a quantum 
critical point |22, H3- Indeed, the T c (x) line in Fig. 
n must terminate at a quantum critical point (QCP), 
which by continuity must be in the universality class of 
the (3 + z)-dimensional XY model, the imaginary time 
t providing the extra z dimensions [z > 1 being the dy- 
namical critical exponent.) Since D — 4 is the upper 
critical dimension for XY-type models, our QCP sits ei- 
ther right at (z = 1) or above (z > 1) its upper critical 
dimension, and thus we expect mean field critical behav- 
ior, possibly with unimportant logarithmic corrections if 
z = 1 [22J, |23j. As indicated in Fig. when crossing 
the finite-temperature transition close to the QCP one 
still encounters the classical fluctuation region, but its 
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width is now much reduced, and it is likely invisible in 
experiments. 

There is a simple consistency check for the above sce- 
nario. If the underdoped region is indeed controlled by 
the (3 + z)D-XY point then there exists a simple scaling 
relation between T c (x) and p s (x,0) which reads [22T E3| 

T c ocp z s /{d - 2+z) . (3) 

For d = 2 we recover the Uemura scaling, irrespective 
of the value of z. In d = 3, as it appears to be the 

z I ( 1 I z\ 

case in YBCO, we get T c cx p s ' y , consistent with the 
experimental observation of Eq. if 1 < z < 2. 

It would thus appear that general arguments based on 
the proximity of the underdoped cuprates to a putative 
(3 + z)D-XY QCP naturally explain items (i) and (ii) 
in the foregoing list. The scaling analysis leading to Eq. 
(|3J), however, only holds when the "bare" parameters of 
the theory describing the critical degrees of freedom ex- 
hibit no significant temperature dependence. In cuprates, 
quasiparticles in the vicinity of d-wave nodes give rise to 
a linear temperature dependence in p s which is likely to 
modify the scaling. In order to address this issue and 
item (m) one needs to go beyond the general scaling ar- 
guments and consider a specific model. In the rest of 
the paper, we formulate and study a particular version 
of the quantum XY model. We show that, when nodal 
quasiparticles are included through the effective param- 
eters of this model, it leads to a phenomenology that is 
consistent with the data. 

The simplest model showing the XY-type critical be- 
havior is given by the Hamiltonian 

ij ij 

Here hi and (pi are the number and phase operators rep- 
resenting Cooper pairs on site of a cubic lattice and 
are quantum mechanically conjugate variables, [hi, ipj] = 
iSij. The sites do not necessarily represent individual 
Cu atoms; rather one should think in terms of "coarse 
grained" lattice model valid at long lengthscales where 
microscopic details no longer matter. Classical [2(| 
and quantum 0, 113 versions of the XY model have 
been employed previously to study phase fluctuations in 
the cuprates. 

The first term in Hxy describes interactions between 
Cooper pairs; we take 

Vi^USij + il-Sij)—^— (5) 

1*2 T j\ 

The second term in Hxy represents the Josephson tun- 
neling of pairs between the sites; Jy = J for nearest 
neighbors along the a and b directions, J' along c, and 
otherwise. In the absence of interactions, J clearly must 
be identified as the physical superfluid density. We thus 
take 

J = Jo- aT (6) 



with a = (21ti2/'k)vf /va, as in a BCS d-wave supercon- 
ductor. The T-linear term describes suppression of the 
mean- field superfluid stiffness by nodal excitations. 

On a qualitative level, the physics of the quantum XY 
model (QJ can be understood in terms of the competi- 
tion between the fluctuations in the local phase (pi and 
charge hi . These are constrained by the uncertainty rela- 
tion Aifi ■ Arii > 1 which implies that interactions, which 
tend to localize charge, also necessarily promote phase 
fluctuations, which then erode the superfluid density. Ul- 
timately, for sufficiently strong Vy, a superconductor- 
insulator (SI) transition takes place. Charge becomes lo- 
calized, and an insulating pair Wigner crystal is formed 
[25|. The latter can be viewed as a superconductor with 
completely disordered phase. In the rest of the paper we 
assume that the strength of interactions increases with 
underdoping and use V{j to tune our model across the SI 
transition. 

To gain quantitative insight into the behavior of p s 
in the XY model (QJ), we employ two complementary 
approaches: the self-consistent harmonic approximation 
(SCHA) [2^, , valid for weak interactions, and an ex- 
pansion in the small order parameter |2JJ, |30j , valid for 
strong interactions in the vicinity of the SI transition. 

In the SCHA one replaces Hxy by the "trial" harmonic 
Hamiltonian 

#har = ^ ™ l Vl 3 "J + \ K V ^i-Vjf- (7) 

ij (ij) 

The constants Kij — K (K 1 ) are identified as the 
renormalized a6-plane (c-axis) superfluid densities [2t| . 
and are determined from the requirement that £"har = 
(HxY)har be minimal. The cuprates are characterized by 
small values of the anisotropy ratio 77 = J'/ J. Anisotropy 
essentially interpolates between the cases d = 2 and d = 3 
and thus profoundly affects the approach to the thermal 
phase transition. However, for the quantum phase tran- 
sition the effect is very weak. We illustrate this point 
below by solving the limiting cases rj = 1 (isotropic 3D 
superconductor) and 77 = (decoupled 2D layers) analyt- 
ically and the intermediate case < rj < 1 numerically. 

The trial Hamiltonian i/har is quadratic in hi and ipj 
and can thus be easily diagonalized, 

ffhar = 5Z?iw q (a^ q + -), fiw q = 2y/KZ^V v (8) 
q 

where Vq is a Fourier transform of Vij and Z q = 
Y =1 sin 2 (g AI /2), d = 2,3. For short range interactions 
Vq — > const as q — > 0; we have w q ~ q, i.e. an acoustic 
phase mode. For Coulomb interactions, Vq ~ l/q 2 as 
q —>■ 0; we have w q — * u) p \, i.e. a gapped plasma mode. 
Simple power counting then shows that at low T the con- 
tribution from the phase mode to the superfluid density 
is 

x ph / T d+1 , short range interaction , . 
\ e _Wpl / T , Coulomb interaction 
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In either case the low-T behavior of p s will be dominated 
by the quasiparticle contribution included via Eq. JSJ. 
However, as we demonstrate below quantum fluctuations 
of the phase lead to strong renormalizations of both the 
T = amplitude of p s and the slope a. 

Using (cos(<p, - ¥?j))har = exp - <pj) 2 )hax], an 

identity valid for harmonic Hamiltonians, we obtain 

#har - (#X Y )har - ^KS - Je'^, (10) 

with the parameter yS = (2dN)^ 1 ^/V^Zq describ- 
ing the aggregate strength of interactions. Minimizing 
-Khar with respect to K we find 

K = Je-^ 1 ^ (11a) 
~ J(l-y/S/J), (lib) 

where the last expression approximates the exact solution 
over much of the regime of interest (see inset to Fig. 

The exact solution of Eq. I|lla|) exhibits a first order 
transition at S w 0.541 J. This is an artifact of the 
SCHA; close to the transition phase fluctuations become 
of the order of tt/2 and the harmonic Hamiltonian JSJ is 
no longer a legitimate approximation to Hxy- Below we 
devise an interpolation formula based on Eq. lllbjl that 
will represent an acceptable solution everywhere except 
very close to the SI transition. 

To obtain the leading temperature dependence, we 
substitute J — Jq — aT into Eq. (jirbfl and expand to 
leading order in T: 

Ps (x,T)~J (l-^S^-aT(^-^^S/j^j. (12) 

This is our main result. As expected, both the am- 
plitude and the slope are reduced by quantum fluctua- 
tions. Crucially, we observe that the T — amplitude 
decays faster than the slope. In particular, for \J S/ Jo 
not too large the above expression is consistent with the 
experimentally observed behavior Eq. @ if we identify 
y ~ (1 — \ \f S/ Jq). In the language of Ref. [2l| the pa- 
rameter y can be interpreted as the quasiparticle charge 
renormalization factor. 

If we follow Lee and Wen Q and determine T c as the 
temperature at which the superfluid stiffness vanishes, 
then Eq. (|12fl implies, to leading order, that p s (x,0) ~ 
Tj?, in agreement with the empirical relation Eq. l|T]). 
The agreement with the scaling result (J3J) is, however, 
entirely coincidental since our description of the super- 
fluid density in the SCHA involves an interplay between 
nodal quasiparticles and noncritical quantum phase fluc- 
tuations. 

For arbitrary anisotropy < rj < 1, the SCHA yields 
a pair of equations for K and K' with structure similar 
to Eq. (| 1 1 a|) . These are easily solved numerically |3lj . 
and we give some representative results in the inset to 
Fig. 12)3. Inspection of this data reveals that all the char- 
acteristic features of the d = 2, 3 limiting cases remain 
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FIG. 2: Panel (a) shows SCHA result, the critical theory Eq. 
114H . and the interpolation discussed in the text (dashed line) 
for the case d = 3 with short-ranged interactions. Inset: solid 
line is the exact solution of Eq. (Illal and the dashed line 
is its approximation Eq. lllbt . Panel (b) displays p s (x,T) 
from Eq. lllbt for several values of S. Inset illustrates the 
effect of anisotropy: the curves (bottom to top) correspond to 
n = 0.0, 0.3, 0.6, 1.0. Solid lines are solutions of Eq. IllaH for 
d = 2, 3 and short range interaction (we find S/U « 0.040 for 
d = 3 and S/U ~ 0.058 for d = 2). Dashed lines are numerical 
solutions of the corresponding anisotropic equations. 



in place for general anisotropy. For realistic anisotropics 
7] = 10~ 2 — 10 -3 the results become practically indistin- 
guishable from the d — 2 case. 

In the regime of strong fluctuations it is useful to con- 
sider the grand canonical partition function for Hxy ex- 
pressed in the path-integral representation as a trace over 
bose field <Pi(r), Z = J Dip exp(— 5 /h), with the action 

1 f p 

2 Jo ij 

Following Refs.lHandliJ, we introduce an auxiliary com- 
plex field V'i( T ) to decouple, via the familiar Hubbard- 
Stratonovich transformation, the cosine term in the 
above action. For short range interaction the decoupled 
action is local in the ipi(r) field and the T>tp functional 
integral can be performed exactly, to any order in powers 
of ip and its derivatives. The field i/j assumes the role of 
the order parameter of the SI transition. Keeping only 
terms up to \ip\ and replacing the spatial lattice by the 
continuum yields the desired field-theoretic representa- 
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tion Z = J T>ipexp(—S c ff/h), where 

ScS = J drd d x jr#| 2 + + i|V^| 2 + T^Wf} 

and r = (d/a§)(l - AdJ/U), u = (7da$- 4 /8)J 2 (4d/U) 3 , 
and c 2 = (4da^/3)(C//4d) 3 / J, with a the lattice spacing. 

The above action S s predicts a second order SI tran- 
sition when r changes sign, i.e. when U = U c = Adj. 
Anisotropy again interpolates smoothly between the lim- 
iting cases [31J and gives U c — 4(2 + rj) J. In (3 + 1) di- 
mensions we expect mean field theory to work near this 
transition. In particular, the superfluid density will be 
given by the saddle-point value of the order parameter 
|t/io| 2 = —r/u which yields 

P°( x > T ) = U-§j) (J-U/12) (14) 

with J given by Eq. 10 and doping parametrized by U. 
The Coulomb interaction can also be incorporated in S c g 
by introducing a gauge field, but the analysis near the 
critical point of the resulting action is more involved [sH ] 
and beyond the scope of this Letter. 

The main panel of Fig. EJa) combines Eq. I|14f) with 
the SCHA result adapted to the case of a short range 
interaction, for which S = 0.48 (77/12). The actual 
solution must interpolate smoothly between SCHA at 
small y/S/J and critical theory near the transition. The 
dashed line represents an empirical extension of Eq. I)llb|) 
to K = J[l - yfS/J - X(S/J)\ with A = 0.625, which 
we expect to be very close to the exact solution, as can 
be verified by quantum Monte Carlo or a similar tech- 
nique. This interpolation still exhibits the leading be- 



havior of Eq. (|12fl . consistent with experimental data 
[Hill 13 El as summarized by Eq. {§. 

Our results thus lend further support to the picture 
of underdoped cuprates as superconductors with a large 
pairing gap scale but superfluid stiffness that is severely 
suppressed by Mott physics. In our approach, the lat- 
ter is modeled by the charging energy terms in the XY 
Hamiltonian J3J which significantly renormalize both the 
T = amplitude of the superfluid density p s (x,T) and 
the quasiparticle effective charge reflected by the slope 
of its T-linear term. A key new observation of this 
work is that the systematics of these renormalizations 
matches that found in underdoped cuprates. In particu- 
lar the suppression of the amplitude is faster than that 
of the slope, in agreement with the experimental data 
[Hlil El Hi summarized in Eq. J5J). As illustrated in 
Fig. [2 this behavior persists over a wide range of inter- 
action strengths. Close to the SI transition, the critical 
theory S c s takes over. In this regime p s {x, T) is given by 
Eq. (|14|l . which implies that the slope of the T-linear term 
stops decreasing and in fact begins to increase. Thus, our 
model offers a testable prediction that very close to the SI 
transition the phenomenology of Eq. |J3J) will ultimately 
break down. While the quantitative details of these pre- 
dictions depend on the specifics of our model and the 
approximations employed, the general features are con- 
trolled by the symmetry and dimensionality and should 
be robust. 
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